Edge states generated by spin-orbit coupling at domain walls in magnetic 

semiconductors 
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Electronic states localized at domain walls between ferromagnetically ordered phases in two- 
dimensional electron systems are generated by moderate spin-orbit coupling. The spin carried by 
these states depends on the slope of the magnetic background at the domain wall. The number 
of localized states is determined by a real space topological number, and spin perpendicular to the 
ferromagnetic order accumulates in these localized states at domain walls. These trapped states 
may be observed in experiments that probe either spin density or conduction paths in quantum 
wells. 

PACS numbers: 71.10. Ay, 71.10. Ca, 71.10. Hf 



Spin transport and spin dynamics in semiconductors 
have attracted considerable recent interest. Spin currents 
in response to an applied electric field have been found 
in both hole-doped semiconductors [l| and the electron- 
doped two-dimensional quantum well with structure in- 
version asymmetry in the third direction [2j . These spin 
currents are referred to as intrinsic spin currents in the 
sense that they do not depend on relaxation by impu- 
rity scattering, and the resulting "spin Hall effect" pre- 
serves time-reversal symmetry, unlike the ordinary Hall 
effect. However, in semiconductors such as GaAs time- 
reversal symmetry can be easily broken if the sample is 
doped with dilute Mn ions, since the Mn ions develop 
ferromagnetic order even at rather low doping 0, Of- 
The ferromagnetic order couples to itinerant electrons or 
holes through both the resulting magnetic field and local 
"Hund's rule" coupling, and hence affects the motion of 
charge carriers; an example is the anomalous Hall effect 
measured in experiments on (Ga,Mn)As [jj. 

This work studies how magnetic domain formation 
modifies conduction electron states in magnetic semicon- 
ductors that also have spin-orbit coupling (recall that 
the existence of spin-orbit coupling alone does not break 
time-reversal). In a sample with spin-orbit coupling, the 
conduction electron magnetization near a domain wall 
(the boundary between oppositely directed domains) has 
a component perpendicular to the external field or the 
local ferromagnetic order, due to one-dimensional elec- 
tronic edge states localized at the domain wall. The pair 
of chiral, spin-polarized electron modes at an edge can 
also lead to an unconventional spin-transport mechanism. 
The mathematical origin of these domain wall states is 
very similar to the origin of trapped edge states in p- 
wave superconductors [6j and their fractional quantum 
Hall analogue 0. 

One possible way to detect the localized domain wall 
states we find is through spatial imagingof spin accumu- 
lation using Kerr rotation microscopy Spin accu- 



mulation changes the direction of polarization of incident 
light. If the incident light is polarized along z direction, 
the bulk spin will not rotate the polarization; only the 
domain wall states can induce a nonzero Kerr angle 6k- 
By measuring the change of polarization of reflected light, 
the accumulation of perpendicular spin can be detected 
at the domain wall. Other possible ways to measure the 
domain wall states are through spin-polarized tunneling 
or imaging of electronic conduction paths |l0| . 

Consider an n-doped narrow-gap semiconductor quan- 
tum well with an asymmetric confining potential in the 
z direction (i.e., perpendicular to the well). In this case, 
the conducting electron band has a Rashba spin-orbit 
coupling pJJl term 



H so = a(k x a y - k y a x ) 



(1) 



In typical GaAs wells, a is about 10" 12 eV-m HI] EH, 
which is quite weak. However, in some materials, for in- 
stance HgTe, the Rashba coupling a is about 100 times 
larger than in GaAs[14j |1It|. and since the Rashba spin- 
orbit coupling is linear instead of quadratic with momen- 
tum, the spin-orbit coupling will become the dominant 
term in the Hamiltonian as long as k is small. Further- 
more, if the sample is doped with Mn ions, then time- 
reversal symmetry could be broken by the long-ranged 
order of local moments of Mn ions. Mn ions are acceptors 
in semiconductor, but the physics discussed here can be 
realized if more donators are also doped other than Mn 
ions. Assume for simplicity that the Mn ions develop 
ferromagnetic order in the z, direction; there will then 
be Hund's-rule coupling between the local Mn ion and 
itinerant electrons, i.e., the local energy of an electron 
depends on whether its spin is parallel or antiparallcl to 
the local moment. We neglect until the end of this pa- 
per orbital magnetic effects, as the exchange energy scale 
that determines the Hund's rule coupling is frequently 
dominant. The full Hamiltonian to quadratic order in k 
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FIG. 1: The spin of localized domain wall states depends 
on the slope of the Zeeman field V. In this example, V is a 
function of x alone, and the edge state is spin-polarized along 
either x or —x. 



is then given by 
A- 2 



H = 



2m* 



+ a(k x a y - k y a x ) + V(x)a z 



(2) 



Let us assume for the moment that the Rashba cou- 
pling is big enough to dominate the quadratic term in 
the Hamiltonian ; the legitimacy of doing this will be 
checked later. The Hamiltonian now reads 



H = a(k x a y - k y a x ) + V(x)a z 



(3) 



If V is a constant, this is exactly the Hamiltonian for 
massive Dirac fermions with mass gap V . Suppose now 
V(x) has a domain wall (a change of sign) at x = 0. 
In order to solve the spectrum of this Hamiltonian, each 
electron state is a spinor with components a(x) and b(x): 



_ JkyV 



a(x) 
b(x) 



(4) 



Since k y is a conserved quantity, the Schrodinger equa- 
tion is given by 



a(-d x b 
a(d x a 



k y b) 

kyO) 



(E-V(x))a 
(E + V(x))b 



(5) 



This equation has two groups of solutions: 

a = —b ~ e" f° ~~^" i<il: , E = ak v 



a = b ~ e^> 



'-dx 



, E = —ak y 



(0) 



If V > 0, which means V(x) is changing from negative to 
positive crossing the domain wall along +x, only the first 
solution is a localized state, and it takes spin o x = —1, 
and the velocity of these states are along +y. On the 
contrary, if V < 0, and V(x) changes from positive to 
negative, the localized state is the second solution above, 
with velocity along —y and spin a x = 1. These eigen- 
states of the Hamiltonian indicates that there is a x ac- 
cumulation at the domain wall, with the sign depending 
on the slope of the local ferromagnetic order (Fig. P). 

Note that © is just a Hamiltonian for massive Dirac 
fermions with variable mass V. That the excitation gap 



closes at the domain wall of V is not so surprising, but 
the explicit solution suggests that these massless domain 
wall states are chiral, with spin depending on the direc- 
tion of velocity of these states. However, it is important 
to note that periodicity of the band structure generates 
an additional set of states with opposite chirality, as dis- 
cussed below. The situation here is analogous to the edge 
states of p±ip superconductors, in which the BdG equa- 
tion breaks into two component equations Q. 

So far the k 2 term has been ignored in all the calcula- 
tions. When is this approximation legitimate? The term 
linear in k will always dominate quadratic terms when all 
the k involved are small. The characteristic k involved 
is determined by the localization length of this localized 
state, which is a/V. Hence the k 2 term is expected to 
be a small perturbation when (V/a) 2 /m is smaller than 
a X V/a, if the value of a in HgTe is taken Q 0] , then V 
should be smaller than 10 meV, which is a typical value 
for magnetic semiconductors. The localization length of 
this state is of order 10 nm. 

We now give a linear stability analysis to justify the 
above heuristic estimate. We have also solved model 
band structures for small systems to verify that trapped 
states indeed exist when predicted by the above intuitive 
criterion. Suppose that the k linear term numerically 
dominates the k 2 term; the stability of the localized so- 
lution for a domain wall state can be checked by solving 
the Schrodinger equation to linear order in 1/m (this lin- 
ear stability is in the same spirit as the WKB method 
in quantum mechanics). For simplicity, V(x) is taken as 
a step function V(x) = Vosgn(x) with Vo > 0. sgn(a;) 
is the sign function of x, sgn(x) = +1 when x > 0, 
and V(x) — —1 when x < 0. We assume at first or- 
der in 1/m that the two components of the solution of 
the Schrodinger equation are modified to be 

(1 + f(x)/m) exp(- / dx V{x)/a) 



ip(x,y) 



= e k vV 



-(1 + g{x)/m) exp(— J dx V(x)/a) 



(7) 



and that the energy is slightly changed to E = ak y + SE. 
As long as / and g have solutions that are polynomials 
of x, the localized states are still stable. After solving to 
first order in 1/m, / and g have the following form: 

f{?) = -o(k y \x\ + hgnix)) + C 



V 1 

g(x) = — (k y \x\ - -sgn(x)) + C 



(8) 



These are solutions with SE = l/m(V 2 /a 2 + k 2 ). These 
solutions suggest that the chiral localized states are sta- 
ble even when there is a numerically small k 2 term. 

The chiral domain wall localized states derived so far 
are not the only states localized at the domain wall, al- 
though they are the most obvious ones. Had the original 
model been defined on the 2d square lattice, the Rashba 
model would read 

H = — t(cosk x + cos k y — 2) + a(smk x a v — s'mk y a x ) 
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The existence of an additional set of states is familiar 
from lattice simulations of chiral fermions: here, the re- 
sult is an additional set of trapped states at a different 
momentum from the original states, as now shown. 

The lattice model coefficient a has different units than 
in the continuum theory. If the continuum spin-orbit cou- 
pling is 10~ 12 eV-m , and the lattice constant is 10~ 10 m, 
the lattice model coefficient a is about lOmeV. The local- 
ized states in previous paragraphs were obtained from lin- 
earizing the Hamiltonian around momentum (0, 0). How- 
ever, the Hamiltonian 10 can always be expanded around 
an arbitrary momentum ko x by defining the slow wave 
function ip(x) from ip(x) — e lko:cX ip(x). 

After linearizing at momentum ko x , the Hamiltonian 
JHJ) becomes 



(9) align spin along +x. These states are only localized when 
s/ a 2 cos 2 ko x — t 2 sin 2 ko x is a real number, which means 



H = tk x sin ko x — t cos k y — i(cos ko x — 2) + Va z 

+a(k x cos ko x a v — sin k y a x ) + a sinko x a v . (10) 



Here {k x , k v ) is the momentum of the slow mode ip. If we 
ignore the last term asmko x a v first (the legitimacy will 
be discussed later), the resulting Schrodinger equations 
for the two-component spinor ip = (a, b) T read 



■ it sin ko x d x a 
= (E + t cos ko x 



— a cos ko x d x b — a sin k y b 



■ t cos k y — 2t — V)a, 



—it sin ko x d x b + a cos ko x d x a — a sin k y a 

= (E + tcosk Qx + t cos k y - 2t + V)b. (11) 

For simplicity, assume that V > (<)0 when x < (>)0 
These two equations have localized solutions as follows: 



\/ a 2 cos 2 k 0x - t 2 sin 2 k 0x . 
E = ; sin k v 

COSfcoz 



-t cos ko x — t cos k y + 2t, 



b = 7a ~ exp( / dx 



V 



V a 2 cos 2 kQ X — t 2 sin 2 k 



O.r 



2Re( 1 ) 2Jm( 7 ) 1 - | 7 | 2 

(cr > = : — lo , (cr y ) = - ; , ,„ , {a } 



i + M 2 



V a 2 cos 2 fc 0a; - t 2 sin 2 k 0x - it sin k 0x 

7 = ; • (12) 

a cos kQ X 

We see for the localized states sin ko x is about a/t, which 
is generally very small, therefore the last term in 11U|) 
is a small perturbation to the whole Hamiltonian. The 
stability of these localized states when k 2 and last term 
in (|10f) is taken into account can be proved in the same 
WKB manner as we derived formula Q. 

As a check, we can assume ko x = and obtain the 
same answer as @. All these localized states tend to 



| tan/cozl < a/t. 



(13) 



However, in order to make sure the approximate solutions 
in (|12l) are approximately orthogonal with each other, 
ko x cannot be taken as continuous value. The inter- 
vals between each ko x should be at least 2irV/a, which 
is the inverse of localization length. Therefore the to- 
tal number of localized channels close to energy Ef = 
is N w 1 + a 2 /(2irtV). There are more than one lo- 
calized channel if a 2 /(tV) > 2tt, which can possibly be 
satisfied in systems with large spin-orbit coupling, like 
HgTe. States with ko x within the intervals {—a/t, a/t) 
and (it — a/t, tt + a/t) satisfy the condition l|13l) . but ac- 
cording to the energy spectrum in (|12l) , the states close to 
kox = 1" have energy higher than states close to (0, 0) by 
2t, and the whole spectrum width of these localized states 
is about 27ra. If the Fermi energy is fixed at Ef = and 
t > ira, only localized states close to ko x — should be 
relevant. These states have the same spin alignment as in 
Figure. ^ but, because the energy E is a periodic func- 
tion on the Brillouin zone, there is always another branch 
of states with velocity opposite to the ones obtained in 
the continuum states. The velocity is v y = dE/dk y , and 
since E is a periodic function of k y , after integrating over 
the whole Brillouin zone the spin current l/2{s x ,v y } van- 
ishes. For example, if t = 0, E = —as'mk y , states close 
to k y — tt have opposite velocity from states close to 
k y = 0. 

Many types of ed ge s tates are related to topology 
in momentum space [lfj. Examples are the quantum 
Hall eff ect 1 1 71. and some quantum spin Hall models 
El III HIM E2- The usual way to bridge the bulk 
momentum-space top olog y and edge modes is through 
Chern-Simons theory [23J, because Chcrn-Simons theory 
is gauge invariant up to a boundary term, which gives rise 
to chiral boson excitations at the boundary. Momentum- 
space topology is usually obtained in insulators, in which 
the whole Brillouin zone is filled; since physical quanti- 
ties have the periodicity of the Brillouin zone, integrals 
over the whole Brillouin zone arc integrals on a compact 
manifold, which can give rise to a nontrivial topological 
number. Here the situation is slightly different: there 
are localized domain wall modes, but the number of edge 
states at a single domain wall is determined by the num- 
ber of bands with linear (Dirac) spectrum near the wall, 
which results when a band is near a minimum or max- 
imum of the band structure neglecting spin-orbit cou- 
pling. Note that the number of domain walls is related 
to the number of zero crossings of the Zeeman field: to 
count this explicitly, a vector N has to be defined. Let 
the z component of N be N z (x) = V(x), and the in-plane 
part of the vector be Ni n = VV. The number of walls is 
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the Skyrmion number 

n=— d 2 xN ■ (d x N x d y N). (14) 

Here N is the unit vector along N. This counting breaks 
down when states from different walls begin to mix. So 
the number of localized states is determined by this real 
space topological number. The spin polarization in Fig. 
^also has nonzero real space Skyrmion number, which is 
equal to ifT^ . 

The domain wall states affect various physical quan- 
tities. The most direct result is the a x magnetization. 
Also, once an electric field is turned on, charge current 
at the domain wall carries spin current (i.e., the charge 
current is spin-polarized) , with spin maximally polarized 
perpendicular to the ferromagnetic order. If the width 
of the sample (or the length of the domain wall) L y is 
shorter than the mean free length of the electrons, the 
modes at ko x = will induce ballistic spin current con- 
ductance similar to the weakly quantized noninteract- 
ing electric conductance in Id systems: for voltage U, 
I s = hwzV , and the total spin current conductance is 
obtained from summing over all the ko x that satisfy (|13J) . 
When a/t < 1, the result is Ne/Air, with N the number 
of localized channels. 

The existence of spin accumulation and spin-polarized 
currents can be understood from symmetry arguments. 
In a 2d plane with asymmetric confining potential in z as 
well as uniform Zeeman field Vz, there are the following 
symmetries: {TP X , TP y , P x P y }. Here P x means x — > —x, 
P y means y — > — y, and T means time-reversal. Because 
T, P x and P y all flip Vz, any combination of two of them 
should be a symmetry transformation of the system. 

Therefore, spin polarization (s x ) and (s y ) are not al- 
lowed by these three symmetries. Instead, spin currents 
jy = l/2{s x ,v y } and j y — l/2{s y ,v x } are allowed by 
symmetry. However, at the domain wall along y, if we 
assume V{x) is an odd function, the symmetries are 
{P x ,TP y }. Now both (s x ) and j x are allowed by sym- 
metry, consistent with the microscopic calculation given 
above. 

For comparison to experiment, let us estimate the rel- 
ative contribution to the electronic density of states from 
the localized states. For an estimate, suppose that the 
Rashba spin-orbit coupling in continuum theory is 10~ 10 
cV-m (as in HgTe), the effective mass is taken as 0.1 m e , 
and if V is about lOmeV, there can be more than one in- 
dependent localized channels at every domain wall. Then 
the spin density (cr x ) per length at the wall is about 10 s H 
cm^ 1 times the number of channels. 

We also find perpendicular spin polarization at domain 
walls when time-reversal symmetry of the itinerant elec- 
trons is broken by an orbital magnetic field rather than 
by a Zeeman field. The Hamiltonian including spin-orbit 



coupling for a 2D electron gas in magnetic field is 

H = a{{k x - eA x )a y - {k y - eA y )a x ) + B(x)g^ B s z (l5) 

The domain wall can be treated as the edge of two quan- 
tum Hall systems, and as is well known, the electron 
states at the edge are already chiral in the absence of 
spin-orbit coupling. The spin direction of these domain 
wall states can be guessed by taking the square of JTSJl: 
any eigenstate of H is an eigenstate of H 2 , although the 
converse is false. (We also ignored the quadratic term in 
the Hamiltonian, since the spin-orbit part will be enough 
for an intuitive discussion.) The square of (|T5)) is 

H 2 = a 2 {k 2 x + {k y - eA y ) 2 ) + 2eBa 2 a z + 2ad x B^i B s x , 

(16) 

where the Landau gauge A x = 0, A y = J Bdx has been 
taken and g = 2. Far away from the domain wall, where 
the magnetic field is varying slowly, A y can be taken 
as Bx. The first two terms are familiar to all, they are 
just the nonrelativistic particles moving in magnetic field. 
The third term results from the slope of the magnetic 
field, and pushes the spin along x, perpendicular to the 
magnetic field. Detailed calculation of the full Hamilto- 
nian gives qualitatively the same result: spin-orbit cou- 
pling leads to confined, maximally spin-polarized edge 
states in the Zeeman case (which is expected to be more 
relevant to experiment), and also tends to spin-polarize 
ordinary quantum Hall edge states in the orbital case, 
perpendicular to the domain wall. 
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